AD-AU4  552  WISCONSIN  UNIV-MADISON  MATHEMATICS  RESEARCH  CENTER  F/6  12/1 

ON  A  CONJECTURE  OF  C.  A.  MICCKELLI  CONCERNING  CUBIC  SPLINE  INTE— ETC(U) 
FEB  62  R  JIA  DAA629-S0-C-0041 

UNCLASSIFIED  MBC-TSR-233* _ _ _ _ NL _ 


one  FILE  COPY 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 

610  Walnut  Street 
Madison,  Wisconsin  53706 

February  1982 


(Received  November  5,  1981) 


Sponsored  by 

U.  S.  Army  Research  Office 
P.  0.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


Approved  for  public  release 
Distribution  unlimited 


dtic 

electej 


UM  188® 

E 


82  05  18  037 


\ 


UNIVERSITY  OF  WISCONSIN  -  MADISON 
MATHEMATICS  RESEARCH  CENTER 


ON  A  CONJECTURE  OF  C.  A*  MICCHELLI  CONCERNING  CUBIC  SPLINE 
INTERPOLATION  AT  A  BI INFINITE  KNOT  SEQUENCE 

Rong-qing  Jia 

Technical  Summary  Report  #2334 
February  1982 


ABSTRACT 

It  is  shown  that  if  the  knot  sequence  t: 
( i )  For  some 


u  J.A.  uic  rwuuu 

■  «  ['•  T-)  • 


'V- 


satisfies 


-1 
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<  JLirn  inf 
r  i 


i+r+1 

'\+T 


ti+r 


-  t. 


<  lim  sup 
r  i 


S.+r+l 


-  t. 
i+r 


'i+1 


-  t. 


<  m 


and 

fci+1  *  fci 

(ii)  m  :  =  sup  — - —  <  •  , 

•  fcj+1  ” 

then  for  any  given  bounded  sequence  y  6  m(8)  there  exists  exactly  one  cubic 
spline  s  with  knots  t^  such  that 

**  t  for  all  i  6  Z  • 
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SIGNIFICANCE  AND  EXPLANATION 


Cubic  spline  interpolation  provides  a  good  and  handy  method  to 
approximate  a  given  function  or  to  fit  a  given  set  of  points.  However,  such 
an  interpolation  process  does  not  always  converge.  It  is  known  that  the  local 
mesh  ratio  (that  of  the  lengths  of  two  consecutive  intervals)  is  less  than 
the  interpolation  process  works  for  any  given  bounded  data. 

This  paper  continues  such  investigation.  It  is  shown  that  the  above 

V 

restriction  on  the  knots  may  be  relaxed.  Thus,  for  a  wider  class  of  knot 
sequences,  the  cubic  spline  interpolation  can  be  still  applied.  Hopefully, 


this  would  make  such  interpolation  process  more  feasible  in  practice. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


ON  A  CONJECTURE  OF  C.  A.  MICCHELLI  CONCERNING  CUBIC  SPLINE 
INTERPOLATION  AT  A  BIINFINITE  KNOT  SEQUENCE 


set 


Rong-qing  Jia 

m 

1.  lij£?E22uction.  Let  t:  -  (t^)^  be  a  biinfinite,  strictly  increasing  sequence. 


t  :=  lim  t  , 
i*±m 


let  k  be  an  integer,  k  >  2,  ana  denote  by  2  the  collection  of  spline  functions  of 

K  f  t 

3C 

degreee  <  k  with  knot  sequence  t.  Explicitly,  S  consists  of  exactly  those  k  -  2 

k ,  t 

times  continuously  differentiable  functions  on 

1  :=  « 

which,  on  each  interval  coincide  with  some  polynomial  of  decree  <  k.  Let 

**k,t  !=  K.t  n  n,(I)  ' 

i.e,,  the  norined  linear  space  of  splines  for  which 

Is  1^  :»  sup  | s ( t ) | 
tei 

is  finite.  We  are  interested  in  the 

Bounded  Interpolation  Problem  (B,I,P«).  To  construct,  for  given  y  6  m(X),  some 

s  £  for  which 

k,t 

( 1 )  s  I 

£  *  y  • 

We  will  say  that  the  b.I,P.  is  correct  (for  the  given  knot  sequence  t)  if  it  has  exactly 
one  solution  for  each  y  £  m(X). 

In  case  k  -  4  (cubic  spline  interpolation),  de  boor  [2]  showed  that  if  the  local 
mesh  ratio 

At 

m  :*  sup  -7—  (with  At  t  ,  -  t.  ) 

s  At,  •  11+11 

i+/b 

is  less  than  — - — ,  then  the  h.I.P.  is  correct,  A  similar  result  was  also  obtained 
independently  by  Zmatrakov  (71.  The  iMsic  idea  of  (2]  was  the  exponential  decay  law,  which 


could  be  traced  back  to  Cl]*  This  idea  was  developed  in  de  Boor  (3)  and  Micchelli  [6J« 

For  cubic  spline  interpolation,  Micchelli  further  raised  the  following  conjecture  (see  (6); 
p*  23b). 


Conjecture*  If 


(2)  m  <  lim  inf  (At^+  /At^  )  4  lim  sup  (At^+r/At^_j )  <  ro 

r  i  r  r  i 

for  some  m  [1#  ) ,  then  the  B.I.P.  is  correct  for  k  -  4. 

In  hindsight,  it  is  easy  to  see  that  this  conjecture  is  faulty.  This  can  be  clarified 
by  the  following 

Example.  Write  h^  :*  t^+j  -  t^.  Let 


fco  -  0 


(I>n'  h2n+1 


(■j) 211+1  for  integer  n  >  0 


h  „  .  *  h.  for  integer  i  >  0 
-l-i  i 


h  h.  — 

*r  <  11m  inf  (r^— )r  <  lim  sup  (.1-^r  )r  <2  , 

•  -  :  **  i  i  _  i 


r  i  i-1 


r  i  i-1 


*1*1  *  \ 


m  -  sup  - - — -  *  -  • 

-  3+1  j 

Let  ra+  and  m“  have  the  same  meaning  as  in  [4).  Since  the  knot  sequence  t  is 
symmetric  with  respect  to  the  origin,  we  must  have  m+  *  m  •  If  this  B.I.P.  is  correct, 
then  (41  tells  us  that  mfc  <  m,  which  is  a  contradiction. 

The  above  example  suggests  to  us  that  the  condition  m  <  •  should  be  added  to  the 
assumption.  Thus  we  will  prove  the  following 

Theorem  1.  If  a  knot  sequence  t  satisfies  (2)  with  1  4  m  <  (3+/5)/2  and 


then  the  B.I.P.  is  correct  for  k  *  4. 

Remark.  This  theorem  covers  de  Boor's  results  for  cubic  splines  with  bounded  global 

_ L  .  .  ^  :  _ -  -  ^  1 _ 1 _ l_ _ vJ-  »  3+/^  . _ f  .  _ |  f  ->  1  \ 


mesh  ratio  or  with  local  mesh  ratio  < 


(see  [2]  and  (3)). 
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2i2-£22iS-£2E52i2S«£2E.SSSi£-5Eii25,i2K£Ki5Si22*  1,01  «  given  knot  sequence  t. 


(4)  h4  tt  -  Vl,  V  -  hi+1/(hi  ♦  hi+1),  U.  -  V<h.  ♦h^,). 


i  6  Z  . 


It  s  6  »*4>t  satisfies  (1)  for  some  y  6  m(Z),  then 


(5)  \sMt._j)  +  2s1  (tA)  ♦  irs't t.+)>  -  3X. 


yi  ~  yi-i  .  ...  vi+1  '  yi 


+  3Mi  -h— - >  1  6  *  * 

i+1 


Moreover,  it  is  easy  to  check  that 


<X'ti-1,<ti-x) 

six)  -  - 2 -  s'(t  ) 

h  1 

i 


(K-Vi)  (VX)  ..  .  yi.,  ♦  v4 

2 -  S'(t  )  ♦  - - -  + 

h  1  *■ 


*  ~2  (x  ~~^2 - L)l‘ti-*)2  ♦  4(ti-x)(x-ti  |)  ♦  (x-t1-t  )2)^i— Yl-1 

"i  "i 


for  x  8  lti_J,ti) 


2  for  j  -  i  , 


Vi  V  Z,  )  \  for  j  *  i“1  / 


\  for  j  -  1+1  , 

VO  for  j  €  . 

Then  A  is  a  tri-diagonal  Z  x  *-matrix.  For  any  0  «  m(z)  *e  have 

Vi  e  Z,  ♦  2B.  ♦  P.S.+1|  <  X.  181.  +  2IBI,  ♦  wJBI..  -  3161.  , 

showing  lAl  <  3.  Mere,  we  view  A  as  a  mapping  from  m(z)  to  m(Z).  Furthermore, 

VI  6  z,  |X.Bi_i  +  26.  ♦  *4i1+||  >  2|B.|  -  X^BI  -  K.  IBI  -  2|  Bi|  -  IBI.  . 
Hence  ,AB*J>  >  sup{2|B^|  -  IBI,,}  ■  IBI,,.  This  shows  that  A-1  exists  and  IA_1  I  <  i. 

3*  J5S„f£E£2t2£iiim£eoaj£.  The  following  lemma  plays  an  essential  role  in  this  pa 
Lemma  For  any  Knot  sequence  t, 

{tt)  lA_1  ( j»i) I  <  Vi,  j  e  z  , 


A-1  (j,i)A_’  (j+1  ,i)  <  0,  Vi,  j  «  z  . 
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Moreover,  it  t  satisfies  the  followinq  condition:  for  some  integer  r  >  0 

(10)  hi^hj  *  *4  whenever  | i — j  |  >  r  , 

then 

(11)  |A1(j,i)|<(l+  m*1  ♦  A  +  +  on’2 )  ^  whenever  |i-j|  *  r 

Proof*  For  simplicity  we  fix  i  and  write  :*  Since  lA  1  I  <  3, 

t>(  • )  6  m(X)*  By  AA™1  *  1  we  have 


(12) 


A  b .  „  +  2b.  +  li  b.  *  6 .  .  »  / 

3  3+1  ID  \ 


A  . 

3  3-1 


0  for  j  f  l 


We  claim  that 
(13) 

If  not,  then  there  exists  jn  <  i  such  that  |b,  |  >  |b.  -|.  From  (12)  we  have 

w  3q  J0f 


!bjl  <  lbj+1l  for  ail  D  <  i 


lb,  -1*  “  l~(2b i  +  "i  bi  +i,/Ai  I  *  (2lbi  I  '  wi  lb-i  ♦,!>/* 

Jn  Jn  jq  J0  J0  Jn 


0  J0 


Hence 


>  (2  -  P  )|b  |/A  -  |b  |  Ml  +  A  )/X  >  2 1  b .  |  . 

J0  J0  J0  J0  J0  J0  Jo 

Then  by  induction  on  j,  we  can  easily  show  that  |b._j|  >  2|b.|  for  all  j  <  j  * 

j0-j  3  3  0 

1 b  |  >2  |b  |  for  j  <  j  .  This  contradicts  the  fact  b  6  m(X)*  Similar  to  (13),  the 
3  3q  u 

following  also  holds: 

(13')  |b^  |  <  lbj.i  I  £or  *11  3  >  i  • 

Now  (12)  and  (13)  yield  that 


lbj+,l  ■  1 2b.  ♦  X.b.^l/M.  >  (2|bj|  -  XjIb^D/p. 


>  (2-X^  )  |b^  |/Uj  *  |b^|*(1  +  >  2lbjl  for  3  <  i  • 


(14) 

Similarly 

04')  \b.^\  >  2|b.|  for  j  >  i  • 

In  particular,  Ib^l  *  lbjj  and  lbi+1  |  4  |b  |.  In  connection  with  (12)  we  obtain 


(15) 


1  -  AiVl  ♦  2b.  +  p.b.+1  «  I V  b.  ♦  2b.  ♦  p.b.+1|  > 

*  2lt>t  I  -  \lbi_,  I  -  >*ilbi+1l  *  <2  -  \  -  VlbJ  -  |o.| 
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This  proves  (tt)  tor  j  =*  i.  For  j  ?  1,  (8)  comes  from  (14),  04')  and  0  5).  For 
rest  of  the  proof  we  may  assume  j  <  1  without  any  loss*  To  prove  (9)  we  argue 


indirectly*  If  b  b.  ,  >  0  for  some  jn  <  i,  then 

Jo  v1  0 

K  _J  *  l\  b,  _,l  -  l2bi  +  **  b- 

■* "*  Jo  Jo  Jn 


Do  V1 


>  2  b. 


Comparinq  the  above  inequality  with  (13),  we  must  have  bj  *  0  for  all  j  4  jQ.  It 
cause  all  b^  **  0,  which  is  absurd*  Now  we  can  write  down 


Ob) 


-b .  , /b .  *;  2  +  2q!  for  j  <  i 
3  +  1'  3  j  J 


with  q  ^  >  0,  Let  :*  hj+l/hj*  We  deduce  from  (12)  that,  for  j  <  i-1. 


b  2b  +  X  b  X  b. 

.  -Ill  .  Ziil - 2111  .  —  +  1J.J,  .  2  ♦  2q  ♦  q 

bj+1  “j+1bj+1  “j+,  “J+,  bj*,  >+1  >  +  1  ^ 


1  4q»  ♦  3 

2  +  2qj+l(’  2(2+2q! ) )  2  +  2qj+1  4qj  +  4  * 


This  shows  that 
07) 

Let 

OB) 


q3+1 


4q*  +  3 
*j+1  4q^  +  4 


for  j  <  i-1 


pj 


4q!  +4 

—2 - q. 

4q‘  +  3"qj 


It  is  easy  to  verify  that 
09) 

Now  (ib)  and  (19)  give  us 
i-1 

n 

k.j 

It  follows  from  (17)  and  (18)  that 


(20) 


IW 


2  +  2q! 


“  lbk+1/bJ 


J  2 

1  ♦  p.  +ri  +  p.  +  p. 


i-1 

n 

k«j 


n  (2  +  2q£)  *  JI  (1 


i-1 

JI 

k=j 


/  2 
♦  pk +  /y  *  pk +  V 


(  21  ) 


i-1  i-1 

n  p  **  it 

k*3  * 

f4qk  + 

4 

_ f  i  \ 

1-1 
«  n 

k*3 

f4qk 

+4 

w + 

3  qk  * 

~+~  3 

4q*_ 

4q  ? 
1- 

1 ♦  4 

4q^ 

+  3 

i-1 

n 

k*j 

3  i_1 

Kn 

k“3 

4q!  “ 

l- 

1 +  3 

#  *^-1 

+  4 

qk  * 

+  3 


1  qkJ 


3  l 
^k  *  4  *  TT 
D 


the 


would 
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If  C  satisfies  (8)  and  |i-j|  >  r,  then 


(22) 


i-1 


n  p 


>  l  .  'll  >  1 . A  .  ,-li-jl 

A  W  A  •%  "a 


k  4  h,  4  3  *0 


k-j  ~  -  **j 

Therefore  lemma  1  will  be  proved,  once  the  following  lemma  is  established: 

Lemma  2.  Suppose  P<i#«***Pn  and  p  are  nonnegative  real  numbers  with  pn 

pip2***pn*  Then 


n  n  f  p.  *  f\  +  p.  +  p^)  >  (i  ♦  p  +  f\  +  p  +  p2)n  • 
.1  1  1  1 


i-1 


Proof.  Let 


n  /  2 

p(p1 #••• #pn-1 #pn>  «■  n  (i  +  +  /i  +  ♦  pt> 


n 

We  want  to  determine  the  minimnm  of  the  function  F  under  the  constraint  D  p  - 

i-1  1 

where  c  is  a  constant,  c  «  pn.  If  some  p^  >  (2+2p)n,  then 


F<P 


Hence 


n  j - — 

•••»,)  >  (2  ♦  2p)n  >  n  (1+P  +  /1+P  +  p2)n  >  in: 

i-i  nPi 


inf  {F(p  ,***,p  )} 
•  _  •  n 


inf  {F(p1 .•••,pn)} 


inf  (F(p -,***P^))  • 

n  in 

Up  — 


Vi,Pi<t2+2p) 


Thus  there  exists  a  point  (p®#...,p°)  with  II  p°  -  c  such  that  P(p?,...,p°)  « 

■  **  4*1  >  n 

0  0  1-1 

inf  {F(p  #•••,?)}•  To  find  (p  ,e«.,p  )  we  shall  use  the  method  of  Lagrange 
ilpt«c  in  in 


multipliers  and  set 


Th#n  l^-l  o 


•CP1#**#Pn)  s-  F<P1»***Pn)  -  XPi*##Pn 
—  0,  i  —  1,*ee,n;  that  is 


,pi  (p“.—.p") 

n  (i  ♦  p?  +  /if  pj  +  (p^)2 )  •  (i 

di  J  J  J 


2p°  +  1 


iti 

It  follows  that 


2/  4-  P°  ♦  <p°>2 


)  -  X  n  p°  «  0  . 


2p°  ♦  1 


p0(l  ♦_4_=r) 


"  p. 


x  • 


1  +  p 


\  +  / 1  +  ?t  ¥  )2  n  (1  >  f  /  +  f  (?^)2) 
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Therefore 


f(p^  *  t(p^)  tor  i*  *  6  {lf..*,n}  , 


f .  2X  +  1  \ 

x(l  +  ■  ■■  :  r-^r-ar^J 

.  2/1  +  X  ♦  X  1  .  X  -  1 

f(x)  - - 7= - 7  “  2  '/"•  '  i 

1  +  x  +  rl  +  x  +  x  2r1  +  x  +  x 


An  easy  calculation  yields 

(  24 )  f '  ( x )  *  - 3-*  1  ~~2  '  3/2  >  0  f  or  x  >  0  • 

4(1  +  x  ♦  x  )  ' 

This  shows  that  f  is  strictly  increasing  on  (0,“).  Thus  (23)  and  (24)  give 

0  0 
P1  *  *  pn  =  p  * 

This  ends  the  proof  of  lemma  2,  Also  the  proof  of  lemma  1  is  complete. 


4.  The  proof  of  theorem  1.  By  the  hypothesis  (2)  there  exist  a  positive  integer  r 

=.= = =»  m  *  * m  s*  amm ss  a  = *=t  ic 

3+/5 

and  a  real  number  mQ  with  m  <  mQ  <  ^  such  that 

h^/h^  <  -j  ^  whenever  |  i— 3  |  >  r  • 

Then  by  lemma  1 , 


(25)  Vi,  j  6  Z,  | A~ 1 ( j , i ) |  < 


2*,i_3 1  if  |i-j|  <  r 


(1  ♦  mQ1  +  /l  +  nig1  +  ,  if  |i-j|  >  r 


M  :=  (m^)  < 


yi  ‘  Yi-1  Yi+1  '  Yi 

c  »  3\.  •  - - - +  3u.  • - ,  i  e  z 

i  X  hi  1  *i+1 


Then  it  follows  that 


(26)  s'(t  )  *  £  A~'(j,i)ct  =  i  A~  (j,i)c.  +  £  A-1 ( j , i )c . 

3  ie*  |i-j|<r  1  )i-jl>r  1 


By  the  hypotheses  of  theorem  1  we  have  the  following  estimates  for 


•7- 


3 

i 


(27)< 


:  I  <  0  •  lyl  •  (i+M)  •—  <  6M ( 1 +M )  •  »y ■  if  |j-i|  <  r  ; 

1  '  a*  tV  nj 

:  I  <  b  •  lyl  •  ( 1 +M )  •  —  <  bO+M)  •  ly  I  •  if  |j-i|  *  t 

i  1  •  n  n .  u 


write 


/  -i  r  -i  -2  \-i 
0  :=  «  It  +  n>0  +  »1  +  i"0  +  mQ  J 


3+/s 

Then  6  <  1  as  long  as  <  — - — .  Applying  (25)  and  (27)  to  (26),  we  obtain 


'  (t.  )j  <  6MO+M)  •  lyl„  •  1“  •  I  2-|l“j|  +  6(1+M)lyl—  •  *■—  I 


3  I  j_i  I  <r 


j  U'j  I  ** 


<28) 


<  const  •  lyl  •  —  • 

«*  h . 

3 


Furthermore  (6)  tells  us 

max  | s ( x ) |  <  const(h  |s'(t  )|  +  h  •  |s‘(t  )|  +  ■y|«)  • 

t.  <x<t.  J  J 

3-1  3 

which  in  connection  with  (2b)  yields  the  desired  result 

tsl„  <  const  •  lyl^  • 

Our  proof  is  complete. 
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